Phonon cooling of nanomechanical beams with tunnel junctions 
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We demonstrate electronic cooling of ID phonon modes in suspended nanowires for the first time, 
using normal metal-insulator-superconductor (NIS) tunnel junctions. Simultaneous cooling of both 
electrons and phonons to a common temperature was achieved. In comparison with non-suspended 
devices, better cooling performance is achieved in the whole operating range of bath temperatures 
between 0.1-0.7 K. The observed low-temperature thermal transport characteristics are consistent 
with scattering of ballistic phonons at the nanowire-bulk contact as being the mechanism limiting 
thermal transport. At the lowest bath temperature of the experiment ~ 100 mK, both phonons and 
electrons in the beam were cooled down to 42 mK, which is below the refrigerator base temperature. 



The ability to cool both the electrons and phonons of 
a mesoscopic device below the temperature of its sur- 
rounding bath is potentially valuable for all areas of low- 
temperature physics, and could prove especially useful 
for applications such as nanobolometry in sub-millimeter 
telescopes [TJ H] , quantum computing [3J and cooling of 
nanomechanical oscillators into their quantum mechan- 
ical ground state [3]. In the sub-Kelvin temperature 
range, the thermal coupling between a nanoscale device 
and its thermal bath becomes very weak [TJ |5J [B] , leading 
easily to overheating problems due to even tiny dissi- 
pated power levels fW) from the measurement signals 
and external noise. A cooling method, where the cooling 
power is applied directly to the device can therefore be 
extremely valuable, as it actually takes advantage of the 
weakness of the thermal coupling. 

A tunnel junction can be used for direct cooling at sub- 
Kelvin temperatures, if one of the electrodes is supercon- 
ducting and the other in normal state, i.e. with a normal 
metal-insulator-superconductor (NIS) structure [Zl [8j [9] 
(Fig. 1). If fabricated on bulk substrates, only sizeable 
electron cooling can be achieved due to the strong weak- 
ening of electron-phonon interaction with temperature 
[SJ [ini [TT]. To be able to cool the phonons, the phonon 
thermal conductance out of the device needs to be small 
enough to become the bottleneck for the heat flow. Us- 
ing this idea, cooling of phonons in thin but large insu- 
lating membranes has been demonstrated with large area 
tunnel junction coolers [TU [T3J [T3], but direct phononic 
cooling of a nanoscale (ID) device has only been sug- 
gested theoretically [15]. In this work, we demonstrate 
significant NIS junction cooling of both electrons and 
phonons in suspended ID nanowires to a common tem- 
perature below the bath temperature of the refrigerator, 
and present quantitative results on the phonon thermal 
transport processes involved (ballistic phonon transmis- 
sion). In our ID geometry, phonon cooling is approxi- 
mately two orders of magnitude more effective than in 
the 2D case [14] in terms of the cooling power required 
to achieve a given temperature reduction. 

All suspended nanowires have a length of either 10 or 
20 /Ltm, width 200 or 300 nm and thickness 60 nm, and 
were fabricated to form bilayers of evaporated copper (30 



nm) on silicon nitride (30 nm) using e-beam lithogra- 
phy, vacuum evaporation and reactive ion etching (see 
[TB]). The Cu/SiN wire is connected to the substrate 
by four free-standing bridges with a length 5 /zm, and 
width 150 nm (Fig. [l] (c)). The outer bridges are also 
Cu/SiN bilayers of thickness 60 nm, and connect the Cu 
wire (total length 20 or 30 /im) to wider (1 ^m) su- 
perconducting Al electrodes on the bulk substrate via 
two larger area (0.35 (/im)^) NIS-junction coolers. The 
cooler junctions must be located on the bulk substrate to 
avoid serious back-flow of dissipated heat from the super- 
conductor into the nanowire 1 17| . The inner bridges have 
a composition Cu/Al/SiN (thickness 90 nm), connecting 
Al leads with Cu quasiparticle traps to two smaller (0.05 
(/im)^) thermometer NIS-junctions located on the sus- 
pended nanowire. Since the measured phonon thermal 
conductance from the narrow bridges to the bulk sub- 
strate is approximately an order of magnitude smaller 
than the thermal conductance between the electrons and 
the phonons (detailed discussion below), phonon trans- 
mission becomes the thermal bottleneck. Thus, both the 
electrons and the phonons in the suspended wire have a 
common temperature and can be cooled simultaneously 
in our sample geometry, in contrast to a recent report on 
electron cooling in nanowires [18J. 

Measurements were perfomed in a '^He-^He dilution re- 
frigerator with a base temperature of ~ 50 mK with sev- 
eral stages of flltering in the wires[16 . In the experiment, 
the temperature of the nanowire was measured as a func- 
tion of the bias voltage across the cooler junctions [16]. 
Schematic of the measurement setup is shown in Fig. [T] 
(c). Using this simple measurement, we obtain tempera- 
ture response curves that typically look like Figure [l] (d) 
as a function of the cooler bias and refrigerator bath tem- 
perature Thath- As expected from theory |9J, the maxi- 
mum cooling is obtained at cooler voltage V ^ 2A/e, 
where A = 215 ^eV for our samples. Two main observa- 
tions are immediately apparent: (i) At the lowest bath 
temperature 50 mK, where the nanowire has a temper- 
ature 100 mK at zero cooler bias (due to noise power of 
^ 6 fW radiated from higher temperature parts of the cir- 
cuit), the lowest temperature achieved around V ~ 2A/e 
was 42 mK, a reduction of ~ 60 mK. (ii) The cooler 
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FIG. 1: (Color online)(a) Schematic of the operation principle of an NIS tunnel junction cooler. The superconducting gap 
A filters through only hot electrons from the normal metal. Heat is removed from the normal metal and injected into the 
superconductor, which must be efficiently removed by quasiparticle traps (b) Theoretically calculated cooling curves 

vs. bias voltage with different hot quasiparticle removal (trapping) efficiencies. The curves with less cooling correspond to 
worse trapping efficiency (effectively hotter superconductor), (c) Schematic of the measurement setup and a scanning electron 
micrograph of a typical suspended sample. Yellow color indicates Cu (normal metal), blue Al (superconductor). Cooler 
junctions are located at the edge of the substrate, whereas the thermometer junctions are at the center of the nanowire. The 
Cu on top of the Al leads serves as quasiparticle trap for both the cooler and thermometer junctions, (d) Measured temperature 
of a typical nanowire sample as a function of cooler voltage at bath temperatures 50,160,200,300,400 and 600 mK, dashed line 
corresponds to 42 mK. 



clearly still works at bath temperatures as high as 600 
mK. These are the main results of this work, and in the 
following we elaborate on the physics by comparing sam- 
ples of different sizes and results between suspended wires 
and wires on bulk substrates. 

Cooler samples on bulk substrates were fabricated with 
the same metal film and tunnel junction geometry as 
the suspended samples during the same fabrication run 
[16j . Figure |2] (a) shows the measured cooling (ratio 
of measured temperature to the bath temperature) vs. 
bath temperature at optimal cooling bias for a suspended 
(green) and a bulk (black) sample, both with wire lengths 
20 /im and similar junction properties, whereas Fig. [2](b) 
shows it for longer 30 /im wires. Clearly, the suspended 
samples in both cases show better cooling (~ 20 % im- 
provement at 300 mK) extending to much higher bath 
temperatures (up to ~ 600 mK), pointing to a differ- 
ence in the dissipation mechanisms between the bulk and 
suspended samples. A further confirmation of the differ- 
ences can be seen by comparing the TiVcooier) curves [T6] . 
In addition, in Fig. [2] (c) we compare cooling for three 
suspended samples in pairs: Two samples have different 
nanowire sizes but the same tunnel junction properties 



(L =20 fim, w =200 nm (black) and L =30 fim, w = 300 
nm (green), Rt '^S kil), whereas the third sample (red) 
has the same size as the longer wire of the previous pair, 
but more transparent tunnel barriers {Rt '^1.7 kfi). It 
is quite clear that the size of the suspended nanowire has 
no effect on the cooling behavior, whereas the tunnel re- 
sistance Rt has a stronger effect, as expected by theory 
fTE\. This size-independence proves that neither 3D [10] 
nor ID [TS] electron-phonon (e-p) interaction limits the 
heat flow in the suspended samples, unlike in bulk cool- 
ers, where a strong volume dependence is observed [S] 
due to the 3D e-p interaction (clearly noticeable by com- 
paring the bulk results in Figs [2] (a) and (b)). We can 
thus deduce that electrons and phonons in the suspended 
wire are in quasiequilibrium (common temperature) and 
are therefore both cooled simultaneously. Moreover, we 
have observed that there are no temperature gradients 
within the wire [TF] so that electronic diffusion [12] is 
not operational either. This means that heat flow is lim- 
ited by phonon transmission (phonon thermal conduc- 
tance), and must be dominated by the interface between 
the suspended wire and the bulk because of the lack of 
wire length dependence. 
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FIG. 2: (Color online) (a) Temperature of a cooled suspended 
(green) and bulk (black) 20 fim long nanowires normalized 
with the bath temperature T/Tbath vs. Tbath- The cooler is 
biased to the optimal cooling voltage, while the bath tem- 
perature is changed. Both samples have an electron gas vol- 
ume Q, — Q.n^ve? and a tunneling resistances Rt ~ 3 kf2 
for suspended and Rt ~ 4.4 kfi for bulk sample. At the 
low temperature regime, cooling efflciency is limited by the 
junctions, whereas at higher temperatures cooling behavior 
is different due to the different dissipation mechanisms (e-p 
interaction for bulk samples vs. phonon transport for sus- 
pended samples), (b) Same, but for longer 30 ^.m long wires 
with O = 0.36 ^vo? and Rt ~ 3 kfi. (c) Same, comparing the 
two suspended samples in (a) (black) and (b) (green) with a 
third suspended sample (red) with O = 0.36 nve? and Rt=1-7 

kn. 

In quantitative terms, the heat flow between the sub- 
strate and the suspended nanowire is determined by 
the power balance condition Pheat — PcooU where the 
heat flow from the surroundings Pheat equals the cool- 
ing power of the junctions Pcooi = "^QcooU where Qcooi 
is the cooling power of a NIS junction[16 . Regardless of 
the limiting heat transport mechanism (e-p interaction 
or phonon thermal conductance), we can generally write 

2Qcooi = A(T,';,, - + p{2Qeooi + IVeooi), (1) 

where the first term on the right describes the the heating 
from the environment (T„^ is the nanowire temperature) , 
and the second backflow of dissipated heat 2Qcooi+IVcooi 
from the superconductor to the normal metal (0 < /3 < 1) 
due to non-equilibrium recombination phonons and back- 
tunneling [205. Here / is the current though the cooler 
junctions and Vcooi the voltage across them, ^ is a pa- 
rameter describing the strength of the heat flow between 
the nanowire and the bath, and the exponent n depends 
on the heat flow mechanism. In the calculation of Qcooi 
we have taken into account a measured broadening of the 
quasiparticle density of states due to lifetime-effects pT] 
and/or two-electron Andreev processes [22] [16], because 
of its strong influence on cooling efficiency [HJ US]. In 
the case of 3D (ID) electron-phonon mediated heat flow 
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FIG. 3: (Color online) (a) Heating experiment: temperature 
of a suspended nanowire as a function of input DC power. 
Red and blue dashed lines correspond power laws n =2.8 and 
n = 6, respectively. The low temperature saturation is likely 
due to absorbed noise power of ~ 6 fW. A = 4.3 pW/K^'* 
for the low-temperature power law n = 2.8. (b) Tempera- 
ture of a suspended nanowire as a function of cooler voltage, 
experimental data is shown as black solid line. The red and 
blue dashed lines correspond to thermal model of eq. [ijwith 
n =2.8 and n = 6, respectively. For the n = 2.8 model, we 
used the value for A determined from the heating experiment 
(no fitting). Inset is zoom-in to the cooling region of the 
low-temperature data in (b). 



A = Sri, where E is the electron-phonon coupling con- 
stant and the electron gas volume (length for ID), and 
n — S — 6 depending on the level of purity of the metal 
film and the phonon dimensionality [HI UHl [H] . In con- 
trast, if the heat flow is limited by phonon transmission 
at the nanowire-bulk boundary (suspended samples), A 
is not dependent on fi, and n = 2 — 6 depending on 
phonon mode and dimensionality [SJ HSl [3^1 [37] . 

To understand the transport process, we performed a 
complementary heating experiment without cooler tun- 
nel junctions but with the same nanowire width w ~ 300 
nm and thickness t = 60 nm (length was L = 24 /im). In 
that experiment we substituted the cooler tunnel junc- 
tions by direct contact between the normal metal (Cu) 
and a superconductor (Nb). These NS-junctions work 
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as good electrical contacts but, on the other hand, as 
nearly perfect thermal barriers so that Joule heating 
power Pin = IV is dissipated uniformly in the normal 
metal and = MTnw " ^Tat/i)- (details in [TB]). Fig. 
[3] (a) shows the result of such an experiment, where the 
temperature of the nanowirc is plotted as a function of 
the heating power in log- log scale. We notice that the 
data is well described by a transition from a power law 
with n = 2.8 at low temperatures to a power law with 
n ~ 6 at high temperatures. The thermal conductance 
G dP/dT at 200 mK is 0.4 pW/K (at 100 mK extrap- 
olated to 0.12 pW/K), which can be compared with a 
value for calculated ID e-p limited conductance [15 (us- 
ing known Cu parameters) at 200 mK, 5.2 pW/K (at 100 
mK 1.3 pW/K), confirming that phonon transport limits 
the heat flow at low temperatures. 

The measured power laws for suspended nanowires can 
then be used in the thermal model of the coolers, Eq[T| 
Figs. [3] (b) and (c) show cooler data in comparison with 
the two power laws, n — 2.8 (red) and n — 6 (blue). 
Good agreement is achieved with the low bath tempera- 
ture data, if a transition from a n = 2.8 power law into 
a n '--^ 6 is assumed, in agreement with the heating ex- 
periment. We would like to stress that simply fitting the 
cooler data with different power laws directly is nearly 
impossible, as the cooler model is quite insensitive to the 
value of n. At bath temperatures Tbath > 0-4:K, n = 6 
fits fairly well the full range of bias values, but is not con- 
sistent with an c-p limited heat flow, as the parameter A 
is much smaller than what is expected from e-p theory 

To understand this behavior, we note that the phonons 
in the wire have a cross-over from 3D to ID behavior 
when the thermal wavelength of the lowest energy trans- 
verse modes At = hct/{2.8kBT) becomes larger than the 
wire thickness and width, which is estimated to take place 
around T ^ 0.4 K using a value Ct = 4300 m/s for the 
Cu/SiN bilayer. This estimate is not far from the ob- 
served transition temperature seen in Fig. |3](b). Thus, 
we believe that in the low temperature regime T < 0.4 K, 
the nanowire phonons are one-dimensional. In the bal- 



listic ID limit with no scattering at the nanowire-bulk 
contact n = 2 is expected [5]. However, in our sample 
geometry the nanowire-bulk contact is abrupt, leading 
to strong scattering and predicted power laws between 
n = 2.5 - 3.5 for 1D-2D scattering [35], or n = 4 - 6 
for 1D-3D scattering [27]. We thus conclude that our ob- 
served n = 2.8 is consistent with 1D-2D boundary limited 
phonon scattering, which is plausible based on the device 
details [TB] . The value of the measured thermal conduc- 
tance (Fig. [3] (a)) per conduction channel G/16 (four 
legs with four phonon branches) is also consistent with 
boundary scattering, as it can be expressed in units of the 
quantum of thermal conductance [5] Go = Tr^kgT/{3h) 
as G/16 - 0.13Go at 200 mK. 

In conclusion, we have demonstrated electronic cooling 
of phonon modes of a suspended nanowire, where scat- 
tering of phonons at the nanowire-bulk contact limits the 
thermal transport. In comparison with non-suspended 
devices, better cooling performance is achieved in the 
whole operating range of bath temperatures. The mini- 
mum temperature reached in our best device was 42 mK, 
starting from an initial temperature of 100 mK. This 
limit is mostly determined by the superconducting mate- 
rial itself (Al) , and thus the minimum temperature could 
be extended below 10 niK by using an additional tunnel 
junction cooler with a lower superconducting gap |28j . 
The advantages of tunnel junction coolers are: (i) ease 
of integration into a nanoscale system and compatibil- 
ity with existing fabrication processes, (ii) simple oper- 
ation by only DC voltage source, and (iii) the ability to 
cool both electrons and phonons simultaneously. This 
last point is significant for ultrasensitive nanobolometers 
[TJ [2], as both the electron and the phonon tempera- 
tures contribute to their performance (noise equivalent 
power), and also for nanomechanical oscillators, where 
the phonon modes need to be cooled. 
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was supported by the Academy of Finland projects No. 
118665 and 118231. 
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I. DEVICE FABRICATION 

The suspended samples were fabricated on single crys- 
tal (001) Si wafers with double-sided 30 nm thick low- 
stress silicon nitride (SiN) films grown by low-pressure 
chemical vapour deposition (LPCVD) and purchased 
from the Microfabrication Laboratory at UC Berkeley. 
A 500 X SOOyitm square was first patterned and etched 
onto the bottom-side nitride, using photoHthography and 
reactive ion etching in CHF3 gas. The remaining SiN 
then served as an etch mask for a crystallographic potas- 
sium hydroxide (KOH) wet etch through the whole wafer 
thickness, producing a 30 nm thick suspended SiN mem- 
brane of size 50 X 50/im on the front side (bulk micro- 
machining). The metallic wires were then deposited on 
the membrane, using electron beam lithography and two- 
angle shadow mask technique, where the the two metals 
aluminium (Al) and copper (Cu) were e-beam evaporated 
from different angles in an UHV system with a base pres- 
sure 10~® mbar. Al was evaporated first (0.2 nm/s) at 
60° angle with respect to the normal of the substrate, af- 
ter which it was thermally oxidized in 10 mbar for 4 min. 
Then Cu was evaporated with rate 0.15 nm/s from 0° an- 
gle. This produced (after lift-off) the narrow nanowire, 
the smaller Al/AlOx/Cu thermometer junctions (typi- 
cal Rt ~ 30 — 50 kfi) in the middle of the wire, and the 
larger Al/AlOx/Cu cooler junctions (typical Rt ^ 1.5—3 
kfi) at the ends of the wire (Fig. 1(c), main text). Af- 
ter metallization, the suspended structure was released 
by reactive ion etching the SiN in a CHF3 plasma (0.1 
mbar for 100 s), where the metal wires themselves served 
as the etching mask. Note that this process also etches 
down the SiN layer on the bulk substrate area around 
the wider leads. 

The bulk samples were fabricated on the same 30 nm 
LPCVD-nitridizcd Si wafers during the same evaporation 
runs, but without any etching processes. 

Samples used in the heating experiment containing di- 
rect superconductor-normal (SN) junctions had a more 
complex fabrication procedure than the coolers, mainly 
because three different materials were used: Cu as the 
normal metal, Al as the superconductor for NIS ther- 
mometers and Nb as the superconductor for the heating 
probes connecting directly to Cu. The thermometer junc- 
tions were still the same Al-A10a;-Cu SIN junctions as in 
the coolers, and thus the first evaporation and oxidation 
steps of Al were the same as for the cooler samples, (60° 
angle, thermal oxidation in 10 mbar for 4 min). Then 
the sample stage was rotated horizontally by 90° and Cu 
deposited from tilt -angle of 60 °. After that, the sample 
was again rotated horizontally by 45° and deposition of 
30 nm of Nb with rate of 0.5 nm/s followed from the tilt- 



angle of 60° . The structural release process was the same 
as for the suspended coolers. A Schematic view of the re- 
sulting heating sample is shown in Fig. [T] Nb was used 
for the heating junctions because of its high supercon- 
ducting gap, which prevents any heat from leaking into 
the superconductor due to multiple Andreev reflections 
[1] . This means that the dissipated power within the nor- 
mal metal nanowire can be accurately determined by sim- 
ply measuring the I-V characteristics of the SNS struc- 
ture, and calculating the dissipated power as Pdiss = IV- 



SN heating 




\ / BHAI 

Probe NIS junctions [ — ] Cu 

■ Nb 



FIG. 1: A Schematic of the sample with two SN heating junc- 
tions. 



II. OPERATION OF THE COOLER 

The basic principle of cooling of a NIS tunnel junction 
is based on the existence of the superconducting energy 
gap A for single particle electronic excitations (Fig. 1(a), 
main text). An electron from the normal metal cannot 
enter the superconductor, unless it has at least energy A. 
A voltage bias V can supply this energy so that at T = 
current can flow if eV > A. However, at finite temper- 
atures electrons follow the Fermi-Dirac distribution, so 
that even at biases eV < A there are some energetic 
(hot) electrons that can tunnel. As only hot electrons 
escape in that case, the temperature of the remaining 
electrons in the normal metal is lowered, and the mag- 
nitude of bias dependent heat flow (power) due to single 
particle tunneling from the normal metal to the super- 
conductor is [2] 

Q{V) ^ 1^ (E ^ eV)gs{E) 

[f^iE - eV, Tn) - fs{E, Ts)]dE, (2) 

where Rt is the tunneling resistance of the junction, 
gs{E) the density of states (DOS) in the superconductor 
and /at.s the Fermi distributions of normal metal and 
superconductor, respectively. Eq. (2) is valid as long as 
both the normal metal and the superconductor remain in 
quasiequilibrium, i.e. as long as the tunneling rate is not 
much larger than the electron-electron scattering rate [5] • 
The optimal cooling power is achieved with a bias close 
to the edge of the superconducting gap V ^ A/e , and 
by increasing the bias further heating starts to dominate 



7 



(Fig. 1(b), main text). It is also clear that to maximize 
the cooling power, the superconducting electrode tem- 
perature Ts has to be as low as possible, which can be 
ensured by removing the excess injected quasiparticles 
as efficiently as possible, for example, by letting them 
quickly diffuse out into another normal electrode (quasi- 
particle trap) pi. For practical purposes, it is useful to 
connect two NIS junctions in series (SINIS), since the 
cooling power is doubled in that case. 

The cooling curves T{V) (Fig. 1(d), main text) were 
measured by slowly sweeping the voltage bias of the two 
cooler junctions in series (SINIS) Kooier and measuring 
the resulting temperature with the second, smaller pair 
of NIS junctions (SINIS thermometer). This was per- 
formed at different refrigerator bath temperatures Tbath- 
In addition, we also performed measurements where the 
cooler voltage bias V^ooior was kept constant at optimal 
cooling point while changing Tbath slowly and measuring 
the temperature response (Fig. 2, main text). This way, 
the performance as a function of Tbath can be determined 
more accurately. Vcooior and Vthorm were both measured 
with a high input impedance differential voltage pream- 
plifier (Ithaco 1201). Measurements of the I{V) char- 
acteristics of the cooler junctions also show consistent 
cooling behavior in comparison with the separate tem- 
perature measurement. 

Because of the extreme sensitivity of the samples even 
to tiny heating powers ^ 10 fW, extra care was taken 
to filter out unwanted external RF noise from the exper- 
iment. The setup uses shielded coax wiring from room 
temperature to 4K flange, where a i?C-filter stage is lo- 
cated. From 4K to 50 mK, the wires are Thermocoax 
cables of length ~1.5 m, with its known good attenua- 
tion properties at high frequencies. In addition, another 
set of RC filters are located at the 50 mK sample stage, 
inside a copper box thermalized to Tbath- 



III. COMPARISON OF THE COOLING CURVES 
OF BULK AND SUSPENDED COOLERS 

The main text discussed the differences between bulk 
and suspended samples in terms of measurements, where 
the cooler bias was kept at optimum value, and bath tem- 
perature was varied (Fig. 2, main text). A further con- 
firmation of the difference between bulk and suspended 
samples is shown in Fig. [2] where measured cooling 
curves (T vs. V) for the bulk (dashes, black line) and 
the suspended (solid, red line) sample of similar values 
of Rt are plotted for three different bath temperatures, 
and compared with the thermal model (Eq. 1 of main 
text) limited either by 3D or ID electron-phonon (e-p) 
interaction in the Cu film (open circles and diamonds, 
respectively) using typical values n — 5 {n = 3 for ID) 
and S = 2.1 • 10^ WK-'^m-^fS]. In addition, for the ID 
model we used c; — 4900 m/s for the longitudinal speed 
of sound in Cu. The only fitting parameter used was /3, 
the fraction of dissipated heat flowing back to the nor- 
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FIG. 2: Cooler sweeps (T vs. V) for a bulk (black, dashed) 
and a suspended (red, solid) sample with Rt ~ 3kfl. Open 
circles were calculated from the thermal model (Eq. 1 main 
text) using typical 3D electron-phonon coupling parameters 
for Cu n=5 and E = 2.1 ■ 10^ WK'^m-^ The fitted value 
for /3 — 0.03 Diamonds are the ID case, with no additional 
adjustable parameters. 



mal metal, with a value /3 = 0.03 kept constant for the 
three bath temperatures and for both models. As we see, 
the 3D e-p model reproduces the bulk data very well for 
all bath temperatures, but cannot explain the suspended 
sample data at all. We stress that since the electron gas 
volumes of the two samples are equal, there are no fit- 
ting parameters left, as we expect f3 to be approximately 
the same for the bulk and suspended samples (both sam- 
ples have the coolers on the bulk substrate). For the 
low temperature range, the suspended nanowire phonons 
are expected to be one-dimensional, thus we should also 
compare the suspended data with the one-dimensional 
e-p model [11 (diamonds. Fig. [2|. It is clear that it can- 
not model the suspended data either (and adjusting /3 
does not improve the fit). The T{V) coohng curves thus 
confirm the conclusion that e-p interaction is not the lim- 
iting dissipation mechanism in the suspended nanowires, 
unlike in the bulk samples. 



IV. DETAILS ON SINIS THERMOMETRY 

In addition to using NIS junction as coolers, they 
were also used as sensitive thermometers because of their 
highly non-linear and temperature dependent current- 
voltage (I-V) characteristics at sub-Kelvin temperatures 
[2115]. In practice, the measurement is usually performed 
by connecting two junctions in series (SINIS), running 
a constant current through the junctions and measur- 
ing the temperature dependent voltage response. When 
constant current biased, the SINIS-thermometer voltage 
^herm is a Sensitive function of temperature only, and 
this dependence can be calculated from the BCS-theory 
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once the tunneling resistance of the junctions Rt and the 
superconducting gap A are known. Since those parame- 
ters were always determined in a separate measurement 
of the I~V characteristics of the junctions, no free pa- 
rameters are left, and the measured SINIS voltage can 
be unambiguously converted to temperature. We al- 
ways checked this to be true by a calibration measure- 
ment, where the SINIS temperature was compared with 
the temperature given by a calibrated RuO thermome- 
ter while the refrigerator temperature was varied. The 
thermometer junctions were current biased with a bias 
resistor R—IQ Gf2 to ensure proper current bias even in 
the subgap, where junction resistance was typically ~10 
MO at low temperatures. In the experiment, two differ- 
ent constant bias current values (/ ^ 10 pA) and (/ '--^lOO 
pA) were used, the low one for low-temperature regime 
and higher one for high-temperature regime. Two values 
were used because the SINIS temperature-to-voltage re- 
sponsivity dV/dT is a strong function of both / and T 
in such a way that the low bias (high bias) value gives a 
better responsivity at T < 0.4 K (T > 0.4 K). In ad- 
dition, the lower bias value was always chosen higher 
than the measured two-electron/lifetime-broadened ex- 
cess sub-gap current, to guarantee good response and no 
dependence on the details of the excess current mecha- 
nisms. Vthorm was measured with a high input impedance 
differential voltage preamplifier (Ithaco 1201). In addi- 
tion, while measuing the I-V characteristics, current was 
measured with a current preamplifier (Ithaco 1211). 

It is clearly seen that by choosing the bias current ap- 
propriately (two examples shown as horizontal lines in 
Fig. [3]), the responsivity dV/dT of the thermometer can 
be adjusted to be best suited for the required temper- 
ature range. Fig. |4] shows an example of the bath- 
temperature-to-voltage response (the bath temperature 
was measured with a calibrated RuO thermometer) of 
the SINIS thermometer in Fig. [3] with bias current val- 
ues 10 pA and 100 pA. With the low bias-current (^ 10 
pA) there is a gain in the responsivity at low tempera- 
tures, while with the higher bias-current (~ 100 pA) the 
responsivity is better at higher temperatures T > 0.4 K. 
Best results with thermometry are obtained by repeating 
experiments with few different bias points for different 
temperature ranges. Note from Fig. [3] that at current 
levels / ^ 1 pA temperature sensitivity is lost. 



A. IV characteristics and subgap conductance of 
NIS junctions 

The I-V characteristics can be understood by a sim- 
ple BCS-based tunneling Hamiltonian theory [3J [S], 
which predicts for a pair of identical normal metal- 
superconductor tunnel junctions (SINIS) 



where Rt is the tunneling resistance of a single junc- 
tion, f^(E) is the Fermi-function of the normal metal 
in quasiequilibrium and gs{E) is the quasiparticle den- 
sity of states (DOS) in the superconductor. Note that 
equation ([3| contains only the temperature of the nor- 
mal metal and not that of the superconductor, and is 
valid as long as the quasiparticle distribution in the su- 
perconductor is also in quasiequilibrium (but not neces- 
sarily the same temperature). gs{E) can be determined 
from the BCS theory, which predicts in the weak coupling 
limit that gsiE) = \E\/^/E'^ - A^, when \E\ > A and 
gs{E) — when \E\ < A However, in real materials 
there are processes that create quasiparticle states also 
within the gap \E\ < A. The easiest and most straight- 
forward way to model these is with the so called Dynes- 
parameter, which was initially realized to model life-time 
broadening due to inelastic scattering (electron-electron, 
electron-phonon) fT. This leads into a DOS of the form 



gsiE) 



Re- 



E + iT 



y/{E + iTy^ - 



(4) 



I{V) 



1 

2eRT 



gs{E)[ME~eV/2)~fN{E+eV/2)]dE, 

(3) 



where parameter F describes the finite life-time (F ~ 
fi/r) of quasiparticle states in the superconductor. 

Figure |5] shows a typical measured suspended SINIS 
cooler I-V curve (black, solid line) in a logarithmic scale 
to highlight the sub-gap regime. As the cooler junctions 
are larger and have thus lower Rt, the current is higher 
than in the thermometer junctions (Fig. [3|. This means 
that the sub-gap current is more easily measurable. The 
dashed, red line in Fig. [sjis a theoretical curve based on 
Eq. [3] that takes into account the broadened DOS by the 
Dynes model (dashed, red line), providing a good fit to 
the data with F/A = 2 • 10~^ as the only fitting parame- 
ter. This value is consistent with what has been observed 
before for thin aluminum films [II [HI |9] . On the contrary, 
without the broadening (dashed, blue line) the sub-gap 
I-Vs can not be fitted. Cooling of the junctions (appar- 
ent as dip around ~ 0.4 mV) is taken into account in the 
theoretical curves by using a voltage-dependent temper- 
ature that was obtained from the separate thermometer 
junctions. 

In addition to the finite life-time broadening, the mi- 
croscopic nature of the finite sub-gap current can also 
originate from the higher order (Andreev) tunneling pro- 
cesses |10l lllj , which give a very similar looking result for 
the I-V curves [12]. For the conclusions of this work it is 
irrelevant which microscopic sub-gap conduction process 
actually takes place, as both of them lead to anomalous 
heating at voltages below the gap and degradation of the 
cooler performance. However, with our observed small 
sub-gap current F/A — 2 • 10""* the anomalous heating 
is a small effect and just barely noticeable in the data at 
lowest temperatures (see Fig. 1(d) in the main text). 
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FIG. 3: A measured sub-gap current-voltage char- 
acteristics of a typical SINIS thermometer {2Rt = 
63 kSl) at different bath temperatures in log-linear 
scale. Solid horizontal lines from top to down cor- 
respond to typical bias currents 100 pA and 10 pA, 
respectively. 




Voltage (mV) 

FIG. 5: Sub-gap current voltage characteristics of SINIS 
Cooler junctions at Tbath ~ 60 mK. Experimental data is 
presented by a solid black line. Dashed lines show numeri- 
cal calculations with broadened DOS F/A = 2 • 10"'* (red) 
and the case where life-time broadening is neglected (green). 
Inset shows the numerically calculated differential conduc- 
tance dl/dV of the measured data (black) and theory with 
F/A = 2 ■ 10"* (red). Rt = 3.3kn. 



B. Temperature measurement in the 
low-temperature regime 

As was shown above, Eq. [3] gives an accurate descrip- 
tion of the measured I{V) curves, even at sub-gap volt- 
ages. This means that if we measure the I{V) curve of 
the SINIS thermometer, we can determine Rt, A and F 
unambiguously, and use Eq. (1) to convert the measured 
SINIS voltage to temperature without fitting parameters. 




Bath temperature (K) 

FIG. 4: Measured voltage-bath temperature response of the 
SINIS thermometer in Fig. [S] with two different values of 
bias current, 10 pA (black line) and 100 pA (red line) (same 
values as the lines in Fig. |3|. Open circles were calcu- 
lated using the single-particle tunneling Hamiltonian, Eq. 
[3] Bath temperature was measured with a calibrated RuO 
thermometer. 

We also point out that by biasing clearly above the sub- 
gap knee (Fig. [3| F does not influence the results at all 
at the temperature range of interest, so that in practice 
only Rt and A fix the temperature calibration. How- 
ever, if one compares the measured responsivity curves 
VsiNis vs. Tbath with the theory from Eq. [s] (Fig. |4|, 
some deviation is clearly seen at the low temperature 
regime T < lOOmK, looking hke a beginning satura- 
tion of VjjiNis- This saturation is not an intrinsic limit 
for the junctions, since we have measured higher values 
of VsiNis (lower temperatures) with the coolers operat- 
ing. We conclude that the observed saturation is most 
likely caused by extrinsic heating power radiated down 
the leads. This is plausible since a nanoscale sample can 
be overheated very easily due to the weakness of dissi- 
pation (electron-phonon, or phonon scattering). From 
the heating experiment (Fig. 3(a) in the main text) we 
can easily see that for our devices an excess power levels 
of 6 fW is enough to heat up the sample to 100 mK 
from Tbath = 60 mK. This kind of power levels are eas- 
ily caused by Johnson noise power radiated down from 
the 4 K stage filters within the bandwidth of our low- 
temperature (below 4K) filters. However, the most im- 
portant conclusion is that this observed saturation does 
not limit our temperature measurement for the coolers, 
as the theoretical self-calibration using Eq. [3] is still valid 
at T < 100 mK, regardless of the excess noise heating. 



C. Lack of temperature gradients in the nanowire 

In the cooling experiments, we measured the temper- 
ature in the middle of the wire, whereas the coolers are 
actually located at the ends of the wire, a distance of 
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FIG. 6: The cooler I~V curve of the sample (black line), 
whose T{V) cooling curves were shown in Fig. 3, main text. 
Tbath = 60 mK. Theoretical curves (Eq. 1) with T as a pa- 
rameter are also shown, with T — 47.3 mK giving the best fit 
at maximum cooling point around V ~ 0.41 V. 



a measurable effect on thermometry, especially in the 
limit Ec > ksT [13], where Ec = /2C^ is the charg- 
ing energy and Cs = 2C + Cq is the total capacitance of 
the junctions and island. Ec can be experimentally de- 
termined by measuring the tunneling conductance spec- 
trum around zero bias in the weak Coulomb blockade 
limit Ec < ksT, where a dip AG develops, depending 
on Ec as dj 



AG 



Ec 



(5) 



where Gt is the tunneling conductance around V — 
without the dip. The measured charging energy at 4.2 K 
for a typical cooler sample is shown in Fig. [7] From the 
measurement we obtain Ec /ks = 20 mK, which is small 
enough not to affect the SINIS thermometry, i.e. analysis 
can be carried out with the simplest BCS-theory calcula- 
tion, Eq. (1), without charging effects. This conclusion 
has been confirmed by a theoretical calculation of the 
SINIS response including charging effects [T3]. 



~10~15 microns away (Fig. 1(c), main text). One might 
wonder whether any temperature gradients will develop 
within the nanowire, and whether the measured tempera- 
ture is therefore unequal to the temperature at the cooler 
junctions. We investigated this question by comparing 
the measured I-V characteristics of the cooler junctions 
with theoretical I-V curves with temperature as a param- 
eter. Fig. [6] It is clear that the cooler I-V is consistent 
with lowest temperature T = 47 mK for this sample, 
which is the same temperature that was measured at op- 
timal cooling for this sample at the SINIS thermometer 
in the middle of the nanowire (Fig. 4(b) main text). 
We conclude that no thermal gradients develop in the 
nanowire. This is consistent with the conclusion that 
phonon transmission at the nanowire-bulk boundary is 
the limiting dissipation mechanism. 



D. Lack of influence of charging effects 

Coulomb blockade (charging effects) due to the small 
capacitance of sub-micron scale tunnel junctions can have 
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FIG. 7; The measured differential conductance spectrum of 
suspended cooler at 4.2K. Conductance is normalized with 
Gt (red, dashed line). The solid line is a one-parameter fit 
to Eq. [5) giving Ec/kB = 20 mK. 
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